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The spin and the chirahty orderings of the Heisenberg spin glass in two dimensions with the 
nearest-neighbor Gaussian coupling are investigated by equilibrium Monte Carlo simulations. Par- 
ticular attention is paid to the behavior of the spin and the chirality correlation lengths. In order to 
observe the true asymptotic behavior, fairly large system size L 20 (L the linear dimension of the 
system) appears to be necessary. It is found that both the spin and the chirality order only at zero 
temperature. At high temperatures, the chiral correlation length stays shorter than spin correlation 
length, whereas at lower temperatures below the crossover temperature Tx , the chiral correlation 
length exceeds the spin correlation length. The spin and the chirality correlation-length exponents 
are estimated above Tx to be usG = 0.9 ± 0.2 and ucG = 2.1 ± 0.3, respectively. These values are 
close to the previous estimates on the basis of the domain-wall-energy calculation. Discussion is 
given about the asymptotic critical behavior realized below Tx ■ 



INTRODUCTION 

Spin glasses (SGs) are the type of random magnets pos- 
sessing both the ferromagnetic and the antiferromagnetic 
couplings, and are characterized both by frustration and 
randomness. Experimentally, it is now well established 
that typical SG magnets exhibit an equilibrium phase 
transition at a finite temperature and that there exists 
a thermodynamic SG phase. The true nature of the SG 
transition and of the SG ordered phase, however, has not 
fully been understood yet in spite of extensive studies for 
years. 

In numerical studies of SGs, much effort has been de- 
voted to clarify the properties of the so-called Edwards- 
Anderson (EA) model [J. Most of these numerical 
works on the EA model have concentrated on the Ising 
EA model. By contrast, many of real SG magnets are 
Heiscnberg-like rather than Ising-like in the sense that 
the magnetic anisotropy is considerably weaker than the 
isotropic exchange interaction jjj, Q . Presumably, theo- 
retical bias toward the Ising model is partly due to the 
simplicity of the Ising model, but partly also reflects the 
historical situation that the earlier theoretical studies on 
the Heisenberg EA model strongly suggested in common 
that the Heisenberg EA model did not exhibit any finite- 
temperature transition Hi S I3> Hi H > apparently making 
it difficult in explaining the experimental SG transition 
within the isotropic Heisenberg model. 

Meanwhile, a novel possibility was suggested by one 
of the present authors (H.K.) that the 3D Heisenberg 
SG might exhibit an equilibrium phase transition at a 
finite temperature, not in the spin sector as usually en- 
visaged, but in the chirality sector, i.e., might exhibit 
a chiral-glass (CG) transition 0- Chirality is a multi- 
spin variable representing the sense or the handedness of 
local noncoplanar spin structures induced by spin frus- 
tration. In the CG ordered state, chiralities are ordered 
in a spatially random manner while Heisenberg spins re- 
main paramagnetic. Refs. 0, H, S Oi Q claimed 



that the standard SG order associated with the freezing 
of the Heisenberg spin occurred at a temperature lower 
than the CG transition temperature at T = Tgc < TcGj 
quite possibly Tsg = 0. It means that the spin and 
the chirality are decoupled on long length scales (spin- 
chirality decoupling). In fact, based on such a spin- 
chirality decoupling picture, a chirality scenario of the 
SG transition has been advanced, which explains the ex- 
perimentally observed SG transition as essentially chiral- 
ity driven 0, ■ Note that the numerical observation of 
a finite- temperature CG transition in the 3D Heisenberg 
SG |l|lH|ll[ll[l3| is not inconsistent with the earlier 
observations of the absence of the conventional SG order 
at any finite temperature. 

Recently, however, in a series of numerical studies 
on the 3D Heisenberg EA model, Tohoku group criti- 
cized the earlier numerical works, claiming that in the 
3D Heisenberg SG the spin ordered at a finite tempera- 
ture and that the SG transition temperature might co- 
incide with the CG transition temperature, i.e., Tqq — 
TcG > 0, Q| . By calculating the spin and the chi- 
rality correlation lengths, Lee and Young also suggested 
TgQ = TcG > 101 . By contrast, Hukushima and Kawa- 
mura maintained that in 3D the spin and the chirality 
were decoupled on sufficiently long length scales, and that 
Tsg < TcG supporting the earlier numerical results. 
The situation in 3D thus remains controversial. 

Under such circumstances, in order to shed further 
light on the nature of the ordering in 3D, it might be 
useful to study the problem for the general space dimen- 
sionality D. The spin and the chirality orderings of the 
Heisengerg SG in dimensions higher than three, i.e., 4D, 
5D and D — oo (an infinite-ranged mean-field SG model), 
were recently studied in Ref.jrg. In the present paper, 
we wish to study the spin and the chirality orderings of 
the low-dimensional system, i.e., the Heisengerg SG in 
2D. 

There were only a few thoeretical and numerical works 
performed for the 2D Heisenberg EA model. On ana- 
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lytical side, there exists a proof that the standard SG 
long-range order (LRO) does not arise at any finite tem- 
perature 01 . The proof, however, does not tell whether 
the CG LRO could exist or not in 2D. On numerical 
side, the domain-wall-energy calculation suggested that 
both the spin and the chirality ordered only at T = 0, 
but with mutually different correlation-length exponents, 
VCG ~ 2.1 > vsG ^ 1.0 — 1.2 0. If this is really the case, 
the spin and the chirality are decoupled on long length 
scale at the T — transition. We note that a similar 
phenomenon has also been reported for the 2D XY SG 
EHHElIi^. Meanwhile, to the authors' knowledge, 
no Monte Carlo (MC) simulation for the 2D Heisenberg 
SG has been reported so far. 

In the present paper, we wish to fill this gap. We study 
both the SG and the CG orderings of the 2D Heisenberg 
EA model by means of a large-scale equilibrium MC sim- 
ulation. The present paper is organized as follows. In §2, 
we introduce our model and explain some of the details 
of the MC calculation. Various physcial quantities cal- 
culated in our MC simulations, e.g., the SG and the CG 
susceptibilities, the spin and the chiral Binder ratios, the 
SG and the CG correlation functions and the associated 
correlation legnths are introduced in §3. The results of 
our MC simulations are presented and analyzed in §4. 
Section 5 is devoted to brief summary of the results. 

THE MODEL AND THE METHOD 

The model we consider is the isotropic classical Heisen- 
berg model on a 2D square lattice with the nearest- 
neighbor Gaussian coupling. The Hamiltonian is given 
by 

H = - ^ JijSi ■ Sj , (1) 

<v> 

where Si = {Sf , Sf , S^) is a three-component unit vec- 
tor, and the < ij > sum is taken over nearest-neighbor 
pairs on the lattice. The nearest-neighbor coupling Jij 
is assumed to obey the Gaussian distribution with zero 
mean and variance . We perform equilibrium MC sim- 
ulations on the model. The lattices studied are square 
lattices with N ^ sites with L = 10, 16, 20, 30 and 
40 with periodic boundary conditions in all directions. 
Sample average is taken over 384 (L = 10, 16, 20) and 
320 {L = 30, 40) independent bond realizations. Error 
bars of physical quantities are estimated by the sample- 
to-sample statistical fluctuation over the bond realiza- 
tions. 

In order to facilitate efficient thermalization, we 
combine the standard heat-bath method with the 
temperature-exchange technique . Care is taken to be 
sure that the system is fully equilibrated. Equilibration is 
checked by the following procedures. First, we monitor 



the system to travel back and forth many times along 
the temperature axis during the temperature-exchange 
process (typically more than 10 times) between the max- 
imum and minimum temperature points. We check at 
the same time that the relaxation due to the standard 
heat-bath updating is reasonably fast at the highest tem- 
perature, whose relaxation time is of order 10^ Monte 
Carlo steps per spin (MCS). This guarantees that differ- 
ent parts of the phase space are sampled in each "cycle" 
of the temperature-exchange run. Second, we check the 
stability of the results against at least three times longer 
runs for a subset of samples. 

PHYSICAL QUANTITIES 

In this section, we define various physical quantities 
calculated in our simulations below. By considering two 
independent systems ( "replicas" ) described by the same 
Hamiltonian J^l, one can define an overlap variable. The 
overlap of the Heisenberg spin is defined as a tensor vari- 
able Qap between the a and f3 components (a, P=x, y, z) 
of the Heisenberg spin, 

1 ^ 

9"/3-^E^-^l? ' K/3 = ^,2/,^) > (2) 
1=1 

where s''^'^ and sf^ are the i-th Heisenberg spins of the 
replicas 1 and 2, respectively. In our simulations, we 
prepare the two replicas 1 and 2 by running two inde- 
pendent sequences of systems in parallel with different 
spin initial conditions and different sequences of random 
numbers. In terms of these tensor overlaps, the SG order 
parameter is defined by 

= u)] . <£= ^i, , (3) 

a,f3—x,y,z 

and the SG susceptibility by 

X^il-N^^ , (4) 

where (■ • •) represents the thermal average and [• ■ •] the 
average over the bond disorder. The associated spin 
Binder ratio is defined by 

Note that gsG is normalized so that, in the thermody- 
namic limit, it vanishes in the high-temperature phase 
and gives unity in the nondegenerate ordered state. 

The local chirality at the i-th site and in the /x-th direc- 
tion Xif^ may be defined for three neighboring Heisenberg 
spins by the scalar, 

Xt,. ^ S^+e^ ■ [Si X Si-e^) , (6) 
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where (p = x, y) denotes a unit vector along the /i-th 
axis. There are in total 2N local chiral variables. We 
define the mean local amplitude of the chirality, x, by 



1 ^ 



(7) 



i=l )[x=x,i 



Note that the magnitude of x tells us the extent of the 
noncoplanarity of the local spin structures. In particular, 
X vanishes for any coplanar spin configuration. 

As in the case of the Hcisenbcrg spin, one can define 
an overlap of the chiral variable by considering the two 
replicas by 



1 ^ (1) (2) 

i—l iJ,—x,y 



(8) 



where Xil^ ^nd xj^-* represent the chiral variables of the 
replicas 1 and 2, respectively. In terms of this chiral 
overlap q-^, the CG order parameter is defined by 



and the associated CG susceptibility by 



XCG = 2N[{ql)] 



(9) 



(10) 



Unlike the spin variable, the local magnitude of the chi- 
rality is not normalized to unity, and is also temperature 
dependent somewhat. In order to take account of this 
effect, and also to get an appropriate normalization, we 
also consider the reduced CG susceptibility xcG by di- 
viding XCG by appropriate powers of X) 



XCG = 



XCG 



The Binder ratio of the chirality is defined by 



(11) 



(12) 



The two-point SG correlation function is defined by 

i 



i a/3 



where f = (x, y) denotes the two-dimensional positional 

vector between the two Heiscnberg spins. The associ- 
ated spin correlation length ^sG is defined, with C{k) = 
C{kx,ky) being the Fourier transform of C(r), by 



&G 




2sin(f) V 



(14) 



Co = (7(0,0), Cm = ^[C{kmin,G) + C{Q,kmin)], (15) 

where kmin = is the minimum nonzero wavevector 
under the periodic boimdary conditions. 

Concering the chirality correlation, we introduce the 
two-point CG correlation function ^^^{r) between the 
two local chiral variables in the /i-th and in the v-ih. di- 
rections, which are apart by the positional vecotr r, by 



^Cg(^") = ^Y.[{Xi,,Xi+r^.)% 



^ ^ ] [(Xi.uXi,uXi+f,i/Xa+f,i/)]' (1^) 



We then define the following three types of chiral correla- 
tion lengths, CcG' ^cg ^cG' depending on the relative 
directions of the local chiral variables with respect to the 
positional vector r, 



CcG — 



2 sin 



(f)\ 



C. 



-1, 



'CG 



2 sin 



c. 



C 



Xm 



CG 



2 sin 



( L ) \ ^xm 



c. 



xo 



-1, 



(17) 



(18) 



(19) 



where the A; = Fourier components C^q, c| q and C^^, 
are given by 

^xo ^^io = l [C'"' (0, 0) + C^nO, 0)] , (20) 

c+o = l[C''no,o) + cy-{o,o)], (21) 

while the k = kmin components C^^, c\m and C+„ are 
given by 

^1 



C^m = o 0) + C^-(0, kr, 



cl 



1, 



xm ^ 2 {f^mim 0) -|- C^^(0, kn 



— '^[C^y {kmim 0) + C'*'^ (^min 7 0) 
-\-C^^{0, kmin) + C^^(0, kmin)]- 



(22) 
(23) 

(24) 



In these configurations, the directions of the two chiral- 
ities are, (i) both perpendicular to the positional vector 
for ^QQ, (ii) both parallel with the positional vector for 
^QQ, and (iii) one perpendicular to, and one parallel with 
the positional vector for ^qq. 
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MONTE CARLO RESULTS 

In this section, we present our MC results. In Fig.l, 
we show the temperature dependence of (a) the SG sus- 
ceptibihty xsG , and of (b) the reduced CG suceptibihty 
XCG (b), on a log- log plot. Recall that the SG transition 
temperature of the present model is expected to be at 
T = 0. As can be seen from Fig. 1(a), our data of XSG is 
consistent with a power-law divergence at zero tempera- 
ture, XSG ~ T~'^^'~', with an exponent 7sg — 1-5 — 1.8. 
In Fig. 1(a), a straight line with a slope equal to —1.8 is 
drawn, which is the best 7sg value determined from the 
other physical quantities analyzed below. Wc note that 
the size dependence of xsG is normal in the sense that 
XSG gets larger as the system size L is increased, and 
that this tendency is enhanced at temperatures closer to 
TsG - 0. 

By contrast, the size dependence of the CG sucepti- 
bihty is somewhat unusual: In the temperature range 
T/J ^ 0.025, xcG gets smaller as L is increased, con- 
trary to the tendency normally expected for a deverging 
quantity in the critical region. At lower temperatures 
T/J'S 0.025, the size dpendence of xcG is changed into 
the normal one, i.e., it gets larger as L is increased. Al- 
though there is no proof that the CG transition temper- 
ature of the model is Tqq = 0, if one fits the data to a 
power-law divergence of the form, xcG ~ T~"'^^ , assum- 
ing TcG = 0, one gets 7cg - 4.0: See Fig. 1(b) The CG 
susceptibility exponent 7cg is then more than twice the 
SG susceptibility exponent 730 ~ 1.8. 

For any T — transition with a nondegenerate ground 
state, the critical-point-decay exponent r/ should be given 
by 77 = 2 — d (d the space dimensionality) , which leads to 
the scaling relation ^ = di' = 21/ (for d = 2) between the 
susceptibility exponent 7 and the associated correlation- 
length exponent ly. Hence, mutually different susceptibil- 
ity exponents for the spin and for the chirality necessar- 
ily mean mutually different correlation-length exponents 
for the spin and for the chirality, i.e., vsg — 0-9 and 
vqg — 2.0. This apparently indicates the spin-chirality 
decoupling. However, one should be careful that the 
present data of the SG and CG susceptibilities alone are 
not enough to completely exclude the possibility that the 
critical behavior crosses over to a different one at still 
lower temperatures, or even the possibility of the occur- 
rence of a finite-temperature CG transition. 

In order get further information, we show in Fig. 2 (a) 
the Binder ratio of the spin gsG, and (b) of the chiral- 
ity gcG- As can be seen from Fig. 2(a), the spin Binder 
ratio for smaller sizes i < 20 exhibits a crossing at a 
nonezero temperature, while the crossing temperatures 
tend to shift to lower temperatures with increasing L. 
For larger sizes L > 20, by contrast, gsG does not exhibit 
a crossing at any temperature studied, always getting 
smaller with increasing L. Such a behavior is indeed con- 



sistent with the expected T = SG transition. Another 
point to be noticed is that gsG for larger sizes exhibits a 
wavy structure at around 0.02 ^ T/J lS 0.03, suggesting 
the occurrence of some sort of changeover in the ordering 
behavior. 




0.01 0.1 
T 



FIG. 1: Temperature dependence of (a) the spin-glass sus- 
ceptibility, and (b) the reduced chiral-glass susceptibility (b), 
on a log-log plot. 

Similarly to the spin Binder ratio, the chiral Binder 
ratio for smaller sizes L < 20 also exhibits a crossing 
at a nonezero temperature, where the crossing temper- 
autres tend to shift to lower temperatures with increas- 
ing L. For larger sizes L > 20, by contrast, gcG does 
not exhibit a crossing at any temperature studied, al- 
ways getting smaller with increasing L. Such a behav- 
ior strongly suggests that the CG transition occurs only 
at T = 0. A closer look of the data reveals that gcG 
exhibits a shallow negative dip at a size-dependent tem- 
perature Tdip{L), which tends to become shallower with 
increasing L. Although such a negative dip is also ob- 
served in the 3D Heisenberg SG, it tends to become 
deeper in the 3D case in contrast to the present 2D 
case0. As argued in R,ef.[l(l|. if a negative dip per- 
sists in the L 00 limit with Tdip{L = 00) > 0, it means 
the occurrence of a finite-temperature CG transition at 
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T = TcG = Tdip{L = oo). In the present 2D case, how- 
ever, the observed Tdip{L) decreases rapidly with increas- 
ing L, tending to T = in the L — > cxo hmit. Indeed, a fit 
of our data of Tdip{L) to the form, TfUpipo) -\- const. /L, 
yields an estimate Tdip{oo)/J = 0.00 ± 0.01. Thus, we 
conclude again that the CG transition of the present 2D 
model occurs only at T = 0, simultancsouly with the SG 
transition. 





FIG. 2: Temperature dependence of (a) the spin Binder ratio, 
and (b) the chiral Binder ratio. The inset is a magnified view 
of the low-temperature region. 

Next, we show our data of the correlation lengths. 
Fig. 3 exhibits the temperature and size dependence of 
the SG correlation length ^sG- One sees from the fig- 
ure that ^sG exhibits a diverging behavior toward the 
zero-temperature transition point. We note that the size 
dependence of ^sg is normal, always getting larger with 
increasing L. 

Figs. 4(a)-(c) exhibit the temperature and size depen- 
dence of the three distinct types of the CG correlation 

lengths, ^cG' '^CG ^^"^ ^CG' defined above. As can be 
seen from the figures, these three chirality correlation 
lengths exhibit very much similar behaviors. They all 
show a diverging behavior toward the T = transition 
point, accompanied by similar size dependence. While 



they all show normal size dependence at lower tempera- 
tures r/J ^ 0.03, increasing with i, they all show un- 
usual size dependence in the range T/J ^ 0.03, decreas- 
ing with L. 



20 



15 



10 



L=10 ' — ^ 

L=16 — * 

L=20 - « 
L=30 

L=40 - - - 



0.02 0.04 0.06 0.08 0.1 0.12 
T 

FIG. 3: Temperature dependence of the spin-glass correlation 
length. 

In order to compare the SG and the CG correlation 
lengths more quantitatively, we plot in Fig. 5 ^sG together 
with , which is the longest one among the three types 
of CG correlation lengths. At higher temperatures, the 
chiral correlation length is shorter than the SG correla- 
tion length. This simply refiects the fact that the short- 
range correlation of the spins is a prerequisite for the 
onset of any chiral order. As the temperature is de- 
creased, ^cG grows much faster than ^sGj then catches up 
and eventually exceeds ^sq at a size-dependent temper- 
ature Tx {L) . This behavior suggests that the chirality 
orders more strongly than the spin. As shown in Fig. 5, 
the linear fit of the log-log plot of the data yields slopes 
i/SG = 0.9 ± 0.2 for ^sG, and vcG = 2.2 ± 0.3 for ^cG, 
the latter being more than twice the former. These val- 
ues are consistent with our estimates of the SG and CG 
susceptibility exponents given above, 7sg = 2i'sg — 1.8 
and 7CG = 2i'cG — 4.0. 

We note that the crossover temperatures Tx {L) where 
^CG exceeds ^sG tend to decrease with increasing L. 
We get Tx/J ~ 0.025,0.019,0.017,0.012 and 0.012 for 
L = 10, 16, 20, 30 and 40, respectively (for L = 30 and 
40, we have extrapolated the data slightly to lower tem- 
peratures to estimate Tx). In particular, if Tx(L) would 
tend to zero as L — > oo, the spin correlation would dom- 
inate the chiral correlation at any nonzero temperature. 
This seems not to be the case here, however, since the 
estimated Tx for our two largest sizes i = 30 and 40 
turn out to be almost the same, Tx/J ~ 0.012. Hence, 
the indication is that, even in an infinite system, ^cG ex- 
ceeds ^SG at low enough temperatures T < Tx ~ O.OIJ, 
where the chiral correlation dominates the spin correla- 
tion. The nontrivial question still to be addressed is the 
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asymptotic critical behaviors of ^sG and of ^cG realized 
below Tx . Although it is difficult to say something def- 
inite in this regime because of the lack of the data, we 
will discuss this point later. 
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FIG. 4: Temperature dependence of the three types of the 
chiral-glass correlation lengths, Eqs.(17)-(19); (a) ^cc (b) 
^CG' and (c) ^Jq. 

In Fig. 6(a), we show the temperature dependence of 
the dimensionless SG correlation length ^sG / L for several 
sizes. Somewhat unexpectedly, ^sg/^ for smaller sizes 



L = 10, 16 and 20 cross almost at a point T/ J ~ 0.022, 
giving an appearance of a finite-temperature SG tran- 
sition. However, ^sg/^ for larger sizes L = 30 and 40 
changes its behavior: They gradually come down, with- 
out a crossing at T/J = 0.022. Hence, the behavior of 
^sg/L for larger L is eventually consistent with the SG 
transition only at T = 0. However, we note that this be- 
havior becomes evident only for larger sizes L ^ 20, and 
one has to be careful in the interpretation of the £,sg/L 
data for smaller lattices. If one looked at the present 
^Sg/L data for smaller sizes only, say, for L < 20, one 
might easily misconclude that the SG transtiion of the 
model occurred at a nonzero temperature T/J ~ 0.022. 
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FIG. 5: Temperature dependence of the spin-glass correla- 
tion length ^sG as compared with the chiral-glass correlation 
length on a log- log plot. 

In Fig. 6(b), we show the temperature dependence of 
the dimensionless CG correlation length ^cG /L for sev- 
eral sizes. Again, ^sG /L for smaller size exhibits a cross- 
ing at a nonzero temperature, while the crossing temper- 
atures tend to shift down to lower temperatures in this 
case. The observed behavior is consistent with the occur- 
rence of a T = CG transition. We note that very much 
similar behaviors are found for the other types of the 
dimensionless CG correlation length, S^^q/L and L 
(data not shown here). 

Having established that both the SG and the CG tran- 
sitions occur only at T = 0, we now analyze the critical 
properties associated with the T — Q transition on the ba- 
sis of a finite-size scaling. In the analysis below, we use 
the data only at lower temperatures T/J < 0.024 where 
the chirality enters into the critical regime exhibiting the 
normal size dependence, and for moderately large sizes 
L > 20. Since the exponent r] associated with the T = 
transition should be zero, we fix Tgc = 0, Tcg = 0, 
VSG = and tjcg = below. 

We begin with the chirality ordering. Fig. 7(a) exhibits 
a finite-size scaling plot of the CG suceptibility, where 
the best value of i^cG is estimated to be vqq = 2.0 ±0.3. 
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The value of pcg is consistent with the one determined 
from the log-log plot of the CG susceptibility. 

In Fig. 7(b), we show a scaling plot of the chiral Binder 
ratio gcG with uqq = 2.0. The data collapse turns out to 
be reasonably good. In Fig. 7(c), we show a scaling plot 
of the dimensionless chiral correlation length ^cG / L with 
VCG = 2.0. The data collapse turns out to be marginally 
good. Hence, our data of the CG susceptibility, the chiral 
Binder ratio and the chiral correlation length seem all 
consistent with vcg — 2.0 ± 0.3. Taking account of the 
corresponding estimate based on Fig. 5, we finally quote 
j^CG = 2.1 ± 0.3 as our best estimate of i^cG- 
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laplse. The quality of the plot does not improve even 
with other choices of z^sG- In Fig- 8(c), we show a scaling 
plot of the dimensionless SG correlation length ^sg/L 
with vsG = 0.9. Again, the data collapse is poor, and 
does not improve even with other choices of I'sG- 
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FIG. 6: Temperature and size dependence of (a) the dimen- 
sionless spin-glass correlation length ^sg/L, and (b) the di- 
mensionless chral-glass correlation length ^Jq/L. The insets 
are magnified views of the low-temperature region. 

Next, we examine the spin ordering. Fig. 8(a) exhibits 
a finite-size scaling plot of the SG suceptibility, where 
the best value of vsg is estimated to be i^sg = 0.9 ± 0.2. 
The value of i>sG is consistent with the one determined 
above from the log-log plot of the SG susceptibility. In 
Fig. 8(b), we show a scaling plot of the spin Binder ratio 
gsG with i^sQ = 0.9. The data collapse is rather poor 
here, however. The wavy structure of the gsG curve ob- 
served for larger L turns out to deteriorate the data col- 
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FIG. 7: Finite-size scaling plots of (a) the chiral- glass sus- 
ceptibility, (b) the chiral Binder ratio, and (c) the dimension- 
less chiral-glass correlation length. The chiral-glass transition 
temperature is TcG = 0. 

Hence, although the SG susceptibility can be scaled 
by assuming i^sg — 0.9, neither the spin Binder ratio nor 
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the SG correlation length can be scaled by the same fsG, 
nor by assuming any other I'sG- Possible reason for this 
is either of the foUowings. 
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FIG. 8: Finite-size scaling plots of (a) the spin-glass suscep- 
tibility, (b) the spin Binder ratio, and (c) the dimensionless 
spin-glass correlation length. The transition temperature is 
TsG = 0. 

i) The width of the critical regime associated with the 
T = SG transition depends on the physical quantities 
measured, being wide for the SG susceptibility but nar- 
rower for the Binder ratio or the correlation length. In 



this case, usg — 0.9 determined above from xsG is a 
true asymptotic exponent, and at low enough tempera- 
tures and for large enough sizes, both gsc and ^sG should 
evnetually scale with i/gQ ~ 0.9. 

ii) The scaling behavior observed in xsG is somewhat 
accidental, not associated with the true asymptotic one. 
At low enough temperatures where ^cG well exceeds ^sGj 
the scaling behavior of xsG is changed into a different one 
characterized by a different i^sG value, which would also 
scale both gsG and ^sG at low enough temperatures and 
for large enough sizes. Unfortunately, the lack of our 
data in the low temperature regime T lS Tx prevents us 
from giving any sensible estimate of the asympotic vsG 
value. If the latter possibility ii) is really the case, there 
even exists a possibility that the spin-chirality decoupling 
eventually does not occur so that one has uqg — ^CG- If 
one trusts our estimate of uqg given above, uqg — 2.1, 
the absence of the spin-chirality decoupling would then 
mean vsg — 2.1. Although we cannot completely rule 
out such a possibility, we also mention that any sign of 
such an equality i^sG = ^^CG was not seen in the available 
data. 

At present, we cannot tell for sure which of the above 
possibilities, i) or ii), is really the case. Nevertheless, 
in view of the scaling observed in xsGi it seems to us 
that the above scenario ii) is rather unlikely, and the 
scenario i) is more plausible. Then, true asymptotic spin 
and chirality correlation-length exponents would not be 
far from i^sg = 0.9 ± 0.2 and j^cg = 2.1 ± 0.3 quoted 
above, which means the occurrence of the spin-chirality 
decoupling in the 2D Heisenberg SG. We note that these 
exponent values are close to the values reported in Ref.0 
for the same model by means of the T = numerical 
domain-wall-energy calculation, i^sq ~ 1.0 ± 0.15 and 
i/CG = 2.1 ± 0.2. 



SUMMARY 

In summary, we performed large-scale equilibrium MC 
simulations of the 2D Heisenberg spin glass. Particular 
attention was paid to the behavior of the spin and the 
chirality correlation lengths. In order to observe the true 
asymptotic behavior, fairly large system size L ^ 20 ap- 
pears to be necessary. We have established that both 
the spin and the chirality order only at T = 0. The 
model seems to possess the two characteristic tempera- 
tures scales T* and {T* > > 0). Around T*, 
the model enters into the chiral critical regime where 
critical fluctuations of the chirality become significant. 
This shows up as a changeover in the size dependence 
of the chiral correlation length or the CG susceptibility: 
Namely, above T*, ^CG and xcG tend to decrease with L, 
whereas below T* they tend to increase with L. The non- 
critical behavior observed above T* might be understood 
as a mean-field-like behavior, as arg ued in Ref. 17]. The 
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chiral correlation length, which stays shorter than the 
spin correlation length at high temperatures, cathes up 
and exceeds the spin correlation length around the sec- 
ond characteristic temperature, Tx ~ O.OIJ. Our data 
in the temperature regime Tx "S T lS T* and for size 
20 5 i 5 40 yields the estimates of the spin and the 
chirality correlation-length exponents, lysG = 0.9 ± 0.2, 
t^CG = 2.1 ± 0.3, respectively. Unfortunately, we donot 
have much data in the interesting temperature range 
T lS Tx, which prevents us from directly determining the 
truly asymptotic critical properties associated with the 
T = transition. While we speculate that the asymp- 
totic spin and the chirality exponents would not be far 
from the values quoted above, the data in the lower tem- 
perature regime T Tx — O.OIJ and for larger sizes 
L 40 is required to settle the issue. 
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